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Abstract 



We consider a problem of estimation for the telegrapher's process 
on the line, say X{t), driven by a Poisson process with non constant 
rate. It turns out that the finite-dimensional law of the process X{t) 
is a solution to the telegraph equation with non constant coefficients. 
We give the explicit law (Pg) of the process X{t) for a parametric 
class of intensity functions for the Poisson process. We propose an 
estimator for the parameter 9 of and we discuss its properties as a 
first attempt to apply statistics to these models. 

Keywords : telegraph equation, inhomogeneous Poisson process, mini- 
max estimation, random motions. 

MSG: primary 60K99, secondary 62M99; 

1 Introduction 

We consider a class of random processes governed by hyperbolic equations. 
These kind of processes have been proposed in the literature to describe mo- 
tions of particles with finite velocities as opposed to diffusion-type models. 
The first contribution in this area goes back to Goldstein (1951). He consid- 
ered the simplest random evolution on the real line where a particle, placed 
in the origin at time 0, moves with two finite velocities ±c changing its cur- 
rent velocity according to a Poisson process of constant rate A. He found 
that the distribution of the position of the particle x at time t, is a solution 
to the telegraph equation, namely 
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where Utt{x,t) = ^u{x,t) and so forth. This model has received a great 
attention in the last decades. Many generalizations have been studied and the 
probabilistic properties of these models have been presented. In particular, 
the finite dimensional and the first passage time laws have been presented in 
an explicit form in a series of paper (see Orsingher 1985, 1990 and 1995, Foong 
1992, Foong and Kanno, 1994). Different generalizations in the presence of 
a finite set of velocities q and Poisson rates Aj (see e.g. Beghin et al 1999) 
as well as motions in two or more dimensions (Orsingher 1986 and 2000, 
Orsingher and Kolesnik 1996, Kolesnik and Turbin 1991, lacus 1995) have 
also been considered. In general, for these models it is quite hard to find 
explicitly the distribution laws of interest. 

In this paper we present a generalization of the model when the underlying 
Poisson process is not homogeneous that are interesting from a statistical 
point of view. Up to our knowledge, this is the first attempt to apply statistics 
in this field. 

The paper is divided into two parts: the first concerns the probabilistic 
analysis of the model and the latter is devoted to statistical estimation. In 
Section || we introduce the model, we study the so called velocity process 
and we derive the following partial differential equation analogous to (Q) : 

utt{x, t) + 2X{t)ut{x, t) = c^u^xix, t) (**) 

where A(-) is the intensity function of the Poisson process (see e.g. Kutoyants 



1998). Then (in §2^) we give conditions on A(-) under which an explicit 
solution to ( p^ ) can be found. It turns out that it is a parametric family of 
intensity measures A(-) = Xe{-)- In the second part of the paper (Section H) 
we consider the problem of estimation of the parameter 6 for this particular 
class of solutions and we present asymptotically efficient -in the minimax 
sense- estimators for 9. 



2 The model and the velocity process 

Consider a particle placed in x = at time t = 0. It can move leftward or 
rightward with finite velocity c. Changes of direction occur at Poisson times. 
This means that the particle will move at constant speed (either +c or — c) 
among two successive Poisson events. We introduce the velocity process V(t): 

y(t) = \/(o)(-i)^w, t>0, (1) 
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where V{0) is a random variable taking values ±c with probability | and 
independent from the Poisson process N{t) with intensity measure A(-) : 



A{t) = [ \{s)ds, t>0, 
Jo 



where A(-) - the intensity function - belongs to C^(]R). The increments 
N{t) — N{s), < s < t, are then distributed according to the Poisson law 
with parameter A{t) — A(s). 

Our main goal is to determine the law governing the process X{t) that 
represents the position of the particle at time t : 

X{t)^V{0) [\-lf^'Ms. (2) 
Jo 

We start with the analysis of the velocity process V{t), t > 0. At first, we 
note that the two probabilities 

p^''\t)^P{V{t)^c) 
=p(y(i) = -c) 

are solutions of the following system of differential equations 

This can be proved by Taylor expansion of the two functions p[^\t) and 
p[ ^\t). The following conditional laws: 

P(V(t) = c\V(0) = c) and F(V(t) = -c\V(0) = c) , 

characterize the velocity process. Their explicit form are as follows. 

Proposition 2.1. 

r* 1 -L p-2A(t) 
P{V{t) = c\V{Q) = c) = l-J^ A(s)e-2^Wds = , (4) 

rt -. _ -2A(t) 

P^V{t) ^ -c\V{0) ^ c) = / A(s)e-2^(^)ds = ^ (5) 

Jo 2 



(3) 
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Proof. We give only the derivation of (^) because the other follows by symme- 
try. Conditioning on V{0) = c implies that: p^'^^{0) = 1, pi'^\o) = — A((0) = 
-A(0), p(-^)(0) = and ^^''^(O) = A(0). From (D we obtain that 

pf^it)~dr^P^ A{t) 
and by simple integration by parts it emerges that: 

pPit) = exp (-2A(t) + log (^] + log(-A(0)) 



,A(0), 

from which follows (^. □ 
Remark that it is also possible to derive (H) by simply noting that 

(oo N 
[j N(t) = 2k 
k=0 

_ 1 + e-2^W 
~ 2 

and then P{V{t) = -c\V{0) = c) = 1 - P{V{t) = c\V{0) = c). 

We conclude the analysis of the velocity process by giving also the co- 
variance function of the couple iV{t), V{s)), s,t > 0. In fact, it can be easily 
proven that the characteristic function of the couple (y(t),V{s)) is, for all 

E (e-^(-)+^/3VW) = cos(«c) cos(/?c) - e-2(A(*)-^(^)) sin(«c) sin(/?c), 
and the covariance function is then 



E(V(s)V(t)) = __^L^E (e'"^(^)+'"^W) 
oaop 



^2g-2|A(t)-A(s)| 



a=P=0 



2.1 Derivation of the telegraph equation 

In order to analyze the distribution of the position of the particle 

P{x,t) = P{X{t) <x) (6) 
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we introduce the two distribution functions F{x, t) =P {X{t) < x, V(t) = c) 
and B{x,t) = P {X{t) < x,V{t) = -c), so that P{x,t) = F{x,t) + B{x,t) 
and W{x,t) = F{x,t) — B{x,t). The function W{-,-) is usually called the 
"flow function" . Next result gives the analogous to (Q) in the case of nonho- 
mogeneous Poisson process. 

Proposition 2.2. Suppose that F{-, ■) and B{-, ■) are two times differen- 
tiable in x and t, then 

Ft{x,t) = -cF,{x,t) - X{t){F{x,t) - B{x,t)) 
Bt{x,t) = cB,{x,t) + X{t){F{x,t) - B{x,t)) 

moreover P{x,t) is a solution to the following telegraph equation with non 
constant coefficients 

g^uix, t) + 2A(t)— t) = c^g^uix, t) . (8) 

Proof. By Taylor expansion, one gets that ■) and ■) are solutions to 
d^) and rewriting system (|^ in terms of the functions W{-, ■) and P(-, ■) it 
emerges that 

' Pt{x,t) = -cW^{x,t) 
Wt{x, t) = cP^{x, t) - 2X{t)W{x, t) ^ ' 

The conclusion arises by direct substitutions. In fact, by the first of system 
(P) we have 

d d 
Pu{x,t) = g^Pt{x,t) = -c—W,{x,t). 

Furthermore 

WUx,t) = —{cP,ix,t)-2\{t)W{x,t)) 
ox 

= cP,,{x,t)-2X{t)W,{x,t) 

= cP,,{x,t) + -\{t)Ptix,t) 
c 

by using respectively the second and the first equation of system (pi). □ 



2.2 The explicit law of the telegraph process 

In this section we give the explicit form of distribution function for a 
particular class of intensity functions. The idea is to reduce equation (|^) to a 
partial differential equation for which the general solution is available. This 
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is done by imposing conditions on the family of intensity functions A(-). Here 
we give only one type of solution that is interesting from the statistical point 
of view. The result presented in the next theorem is interesting in itself. 

Theorem 2.3. Suppose that the intensity function of the Poisson process 
N{t) m (g) IS 

X{t) = Xe{t) = e tanh(^t), OeR. (10) 

Then, the absolutely continuous component p0{-) of distribution is given 
by 

^ \x\ < ct 



po{x,t,c) = -{ ^^^Hdt) 2Vc^W^ ' i-i--" (11) 
0, otherwise 

Proof. We start without assuming ( [T0| ) and by noting that p{x, t) is a solution 
to (H). Moreover, the process X{t) G {—ct, ct) only if there occur at least one 
Poisson event up to time t, thus 

/ct 
p{x, t)dx = 1 - PiN{t) = 0) = 1 - e"^(*) (12) 
■ct 

that gives one of the conditions to solve (^. We now search for solutions of 
the following type: 

v{x,t) = e^(*)M(x,t) . 
Thus the function f (■, ■) satisfies the following partial differential equation 

vttix, t) - vix, t) (X'it) + \^{t)) = cVx(a;, t) . (13) 



For a generic function A(-) a solution to (|13[) not available. If the intensity 
function A(-) satisfies the following ordinary differential equation 

X'{t) + X\t) = e^, t>0, 9eR, (14) 

then equation (|T3|) becomes 

vtt{x, t) - v{x, t)9'^ = C^Va:x{x, t) (15) 

and the solution to it can be written explicitly^. With the initial condition 
A(0) = the solution to equation (|1^) is condition (p!0| ) and then Ag(t) = 

-'^ Remark that, if one solves the problem A'(t) + A^(i) — equation ( p^ ) reduces to 
the standard one-dimensional wave equation v^x = c^Vxx and its solution can also be 
determined but it is of no statistical interest. 
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ln(cosh(6't)). Following Orsingher (1985), by the change of variable s = 
a/ c^t^ _ transform equation ([TsD into the following standard Bessel's 

equation 

1 



whose general integral is 



(0 



Vss + -Vs- ( - ) v = (16) 



c 



The function /fc(x) is the modified Bessel function of first kind and order k 
while Ko{x) is the second type Bessel function of order with the unpleasant 
property that lim^^o+ Ko{x) = oo, so we put B = 0. In terms of pe^-, ■) the 
solution is of the following form 



|x|<ct. (17) 



From (|T2|) it follows that 

ct 



[ pg{x,t)dx = 1 - e-^o^^^ = 1 ^ 

J cosh 



-ct 



cosh(^t) e"^* + e^* 



and so, there is no constant value K that satisfies ([T7|) . To turn around this 
problem, we note that if Iq [^y/cH"^ — x^) is a solution to (|^) so is its partial 
derivative with respect to t. Hence we search for solutions of the form 



Form the following two equalities (see Orsingher, 1995) 

/, (i^^cH^^x^^ dx = ^- (e^* - e-^*) 



and 

-ct 



£ ^/o (^Vc2t2-x2^ dx = c (e^* + e-^*) - 2c 
condition (|12[) implies that 

1 _ e-'^^W = e-^«W (a(1 - \g{t))^ (e^* - e"^*) + 5c (e^* + e''' - 2 



1 - e-at+Jt thus the 

equahty in the above equation is attained by taking A = and B = The 



Observe that for Xe(t) as in (p!0| ) we have 1 — \e{t) 
equahty in the abi 
solution is finaUy 

dth (fVc2 



Pe{x,t) 



\x\ < ct 



that is non negative and satisfies the conditions required. It also verifies 
equation (H) and this is a boring calculus' exercise. □ 



Note that the non-null component of (P) can be written as 



2 cosh(^t) 



{5{x - ct) + 5{x + ct)) + 



eth {^VcH^-x^) 



x^ 



X G [—ct, ct] 



where S{x) is the Dirac delta. By the result of last theorem we can now state 
for the seek of completeness the following result. 

Proposition 2.4. Suppose that the intensity function of the Pois son process 
N{t) m (D IS 

\{t) = \e{t) = etim\{et), 0gM. 



Then, the distribution function 



IS 



Pe{x,t) 



0, 



1 



2 cosh(et) 
1, 



+ / pe{u,t)du, 



-ct 



X < —ct 
—ct < X < ct 

X > ct 



3 Parameter estimation 

We consider the problem of estimation of the parameter 9 for the model 
discussed in § |2.2| . We suppose that the velocity c is known, if it is not, 
one can easily determine without error its value, by observing two successive 
switches of velocity (or just the first switch) and the time between the two 
occurrences. Two possibile schemes of observations can be considered: i) 
one trajectory is observed up to some time T, letting T — oo, and m) n 
independent and identically distributed observation Xi{t) of the trajectories 
are observed on a fixed time interval [0,T], letting n oo. Instead of 
working on the law of X{t) we do inference on 9 via the Poisson process. 
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We use the method of moments to estimate 6. For scheme i), recall that 
EA^(t) = K0{t) thus, if ttt is the number of observed switches up to time T, 
it suffices to find the solution to 

KeiT) = ln(cosh(eT)) = ttt 

that is 

Good asymptotic properties of this estimator cannot be established without 
assuming more conditions on the Poisson process. We give now the optimal 
solution in terms of the second scheme. 

Suppose now to be able to observe n independent copies Xi{t) of the 
process X{t) up to a fixed horizon T, as to say, wc observe the trajectories 
of n particles that do not interact. Denote by tTj the number of switches in 
each replication. Then 

n 

- E 



n 

1=1 



is a consistent estimator of tt = ^e{T) as n ^ oo, moreover y/n{7{n — tt) 
is asymptotically Gaussian A/'(0,7r). The estimator 7r„ - called empirical 
measure - is also asymptotically efficient in the minimax sense (for these 
and many other results see Kutoyants, 1998). In this particular case, an 
asymptotically efficient estimator for exists. In fact, 

TT = ln(cosh(^T)) and 9 = 9{tt) = ^ arcosh (e'^) 

so the estimator of moments is 9n — ^(7r„). By the so-called 5-method we 
desume the asymptotic behavior of ^(7r„). In fact 



Vn{9{7rn) - 9{7r)) = v^(7r„ - n)9'{n) + o(|7r„ 



TT 



where the last term converges to in probability by consistency. Thus, 
V^{9{7rn) — 9{t:)) has, asymptotically, a centered Gaussian law with variance 
given by 



7 - TT 

7r/^'(7r)2- 



1 + coth(7r) 



or better 



This estimator is also the maximum likelihood estimator for 9 that is asymp- 
totically efficient in this regular problem. To conclude, the estimator 9n is 



9 



consistent for 9 and asymptotically Gaussian, moreover it is asymptotically 
efficient for 6 with asymptotic variance given by the inverse of the Fisher 
information in this problem. 
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